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Petrov type D gravitational fields, generated by a perfect fluid with spatially homogeneous energy 
density and with flow lines which form a non-shearing and non-rotating timelike congruence, are 
re-examined. It turns out that the anisotropic such spacetimes, which comprise the vacuum C- 
metric as a limit case, can have non-zero expansion, contrary to the conclusion in the original 
investigation by Barnes [l||. Apart from the static members, this class consists of cosmological 
models with precisely one symmetry. The general line element is constructed and some important 
properties are discussed. It is also shown that purely electric Petrov type D vacuum spacetimes 
admit shearfree normal timelike congruences everywhere, even in the non-static regions. This result 
incited to deduce intrinsic, easily testable criteria regarding shearfree normality and staticity of 
Petrov type D spacetimes in general, which are added in an appendix. 
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1. INTRODUCTION 



The C-metric is a well-known exact solution of Ein- 
stein's vacuum equation with zero cosmological constant. 
The static region of the corresponding spacetime was first 
described by Weyl [2j. At about the same time Levi- 
Civita [3] constructed its line element in closed form, ar- 
riving at essentially one cubic polynomial with two pa- 
rameters as the metric structure function. The C-metric 
is a Petrov type D solution for which at each space- 
time point both Weyl principal null directions (PNDs) 
are geodesic, non-shearing, non-twisting but diverging; 
it thus belongs to the Robinson- Trautman class of so- 
lutions and was rediscovered as such [4]. The label 'C 
derives from the invariant classification of static degen- 
erate Petrov type D vacuum spacetimes by Ehlers and 
Kundt [5!|. The importance of this solution as summa- 
rized by Kinnersley and Walker Q, is threefold. First, 
the C-metric describes a spacetime with only two inde- 
pendent Killing vector fields (KVFs) which can be fully 
analyzed. Next, it is an 'example of almost everything', 
most notably it describes a radiative, locally asymptot- 
ically flat spacetime, whilst containing a static region. 
The C-metric is contained in the class of boost-rotation- 
symmetric spacetimes 0, 1^, which are the only axially 
symmetric, radiative and asymptotically flat spacetimes 
with two Killing vectors. Finally, the solution has a clear 
physical interpretation as the anisotropic gravitational 
field of two Schwarzschild black holes being uniformly 
accelerated in opposite directions by a cosmic string or 
strut, provided that ma < l/-\/27, where the mass m 
and acceleration a are equivalents of the two essential 
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parameters of Levi-Civita @, Q (see, however, the end of 
12.31 for a comment). 

Generalizations of the C-metric have been widely con- 
sidered. Adding a cosmological constant A is straightfor- 
ward, and we will henceforth refer with 'C-metric' to such 
Einstein spaces. Incorporating electromagnetic charge 
|(7p = 4- 5^ is equally natural and leads to quartic 
structure functions . Recently, the question how to in- 
clude rotation for the holes received a new answer ^10l . lTl| . 
avoiding the NUT-like behavior of the previously consid- 
ered 'spinning C-metric' (T2| . IT3| . All these generaliza- 
tions fit in the well-established class V of Petrov type D 
Einstein-Maxwell solutions with a non-null electromag- 
netic field possessing geodesic and non-shearing null di- 
rections aligned with the PNDs [13, [l^ , which reduces 
for zero electromagnetic field to the subclass I?o of Petrov 
type D Einstein spaces and which contains all well-known 
4D black hole metrics. In fact, all 2?-metrics can be de- 
rived by performing 'limiting contractions' [iG] from the 
most general member, the Plebianski-Demianski line ele- 
ment [l7| , which exhibits two quartic structure functions 
with six essential parameters m, a, jgp. A, NUT pa- 
rameter / [Tsj and angular momentum a. A physically 
comprehensive and simplified treatment can be found in 
[1^], also surveying recent work in this direction. 

In this paper we present a new family of Petrov type 
D, expanding and anisotropic perfect fluid (PF) gener- 
alizations of the C-metric. The direct motivation and 
background for this work is the following. 

According to the Goldberg-Sachs theorem [20] the two 
PNDs of any member of are precisely those null di- 
rections which are geodesic and non-shearing. Such a 
member is purely electric (PE, cf. appendix B) precisely 
when both PNDs, as well as the complex null direc- 
tions orthogonal to them, are non-twisting (non-rotating 
or hypersurface-orthogonal (HO)). This is in particular 
the case for the C-metric. As we will show, it implies 
the existence of an umbilical synchronization (US), i.e., a 
non-shearing and non-rotating unit timelike vector field 
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(tangent to a congruence of observers). The importance 
of USs in cosmology was stressed in [2lj. If a congru- 
ence of observers measuring isotropic radiation admits 
orthogonal hypersurfaces, an US exists. Only small de- 
viations from isotropy are seen in the cosmic microwave 
background, and scalar perturbations of a Friedmann- 
Lemaitre-Robertson- Walker universe preserve the exis- 
tence of an US [22]. In general, spacetimes admitting 
an US have zero magnetic part of the Weyl tensor wrt 
it [23!| and thus are either of Petrov type O, or PE and 
of type D OY I [l^. Conformally flat spacetimes always 
admit USs (see e.g. (6.15) in [16|). Triimper showed that 
algebraically general vacua with an US are stat ic |24| . 
Motivated by this result and by his own work [25| on 
static PFs, Barnes f]] studied PF spacetimes with an 
US tangent to their flow lines. He was able to gener- 
alize Triimper's result to Petrov type / such PFs and 
recovered Stephani's results on conformally flat PF solu- 
tions which are either of generalized Schwarzschild type 
or of generalized Friedmann type (so called Stephani uni- 
verses) {2^. The type D solutions were integrated and 
invariantly partitioned, based on the direction of the gra- 
dient of the energy density relative to the PNDs and the 
flow vector at each point. Class I, characterized by the 
energy density being constant on the hypersurfaces or- 
thogonal to the flow lines and thus the only class contain- 
ing Einstein spaces as limit cases, was further subdivided 
using the gradient of ^'2 (cf. ij 12.21 for details). By solv- 
ing the field equations, Barnes concludes that class ID, 
consisting of the anisotropic class I models, solely con- 
tains non-expanding solutions. Hence, these PF solutions 
would not be viable as a cosmological model. However, 
based on an integrability analysis of class I in the Geroch- 
Held- Penrose (GHP) formalism |23|, we found that this 
conclusion cannot be valid and this led to a detailed rein- 
vestigation. 

In this article we construct the general line element 
of the full ID class, comprising both the known non- 
expanding perfect fluid models and the new expanding 
ones, and discuss some elementary properties. We want 
to stress the following point. The full class represents a 
PF generalization of the C-metric in the sense that the 
C- metric is contained as the Einstein space limit. The 
physical interpretation of this fact is however not estab- 
lished. This would require to exhibit this solution for 
small masses as a perturbation of a known PF solution, 
just as the C-metric interpretation of small accelerating 
black holes has been established in a flat or (anti-)de Sit- 
ter background (6l.[2i-[3H. 

However, the mathematical relation with the C-metric 
is useful. As already deduced in [1|, the PF solution is, 
just as the C-metric, conformally related to the direct 
sum of two 2D metrics. The fact that one part is equal 
for the PF solution and the C-metric is helpful in the 
analysis, e.g. we will show that (a part of) the axis of 
symmetry can readily be identified as a conical singular- 
ity, analogous to the defect of the cosmic string present 
in the C-metric. The non-static spacetimes presented 



are exact perfect fluid solutions with only this symme- 
try, and the analysis appears to be within reach. For the 
expanding ID PF models both the matter density ■w{t) 
and the expansion scalar 9{t) can be arbitrary functions. 
This freedom is displayed explicitly in the metric form, 
and makes the solutions more attractive as a cosmological 
model. 

The paper is organized as follows. In section 2 we 
present the GHP approach to class I. We derive a closed 
set of equations, construct suitable scalar invariants, in- 
terpret the invariant subclassification of [1] and start the 
integration. At the end we provide alternative charac- 
terizations for the Einstein space members and identify 
their static regions and USs. In section 3 we finish the 
construction of the general ID line element in a trans- 
parent way, and correct the calculative error of [1] in the 
original approach. Then we deduce basic properties of 
the ID perfect fluid models. In section 4 we summarize 
the main results and indicate points of further research. 
The work greatly benefited from the use of the GHP for- 
malism, which at the same time elucidates the deviation 
from the C-metric. In appendix A we provide a prag- 
matic survey of this formalism for the non-expert reader. 
In appendix B, finally, we present criteria for deciding 
when a Petrov type D spacetime admits a (rigid) US or 
is static. 

Notation. For spacetimes (Af , gab) we take (+ + + 
-) as the metric signature and use geometrized units 
8ttG — c — 1, where G is the gravitational coupling con- 
stant and c the speed of light. A denotes the cosmological 
constant. We make consistent use of the abstract Latin 
index notation for tensor fields, as advocated in (s^ . 
Round (square) brackets denote (anti-)symmetrization, 
Vabcd is the spacetime alternating pseudo-tensor and 
'^cTab... ('CxTah...) designates the Levi-Civita covariant 
derivative (Lie derivative wrt X"-) of the tensor field 
Tab...- One has 

daf^Vaf, dbYa^V[bYa] 

for the exterior derivative of a scalar field /, resp. one- 
form field Ya, and 

X(/)=X"d„/ 

denotes the Leibniz action of a vector field AT" on /; 
when AT" is the x'-coordinate vector field d^i"" we write 
d^i f or f ,^i , and a prime denotes ordinary derivation for 
functions of one variable, f'{x) = dxf{x). However, we 
use index- free notation in line elements ds^ = gijdx''dx^ . 
The specific GHP notation is introduced in appendix A. 

2. GHP APPROACH TO CLASS I 

2.1. Definition and integrability 

We consider Barnes' class I [l], consisting of spacetimes 
(M, gab) with the following properties: 
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(i) the spacetime admits a unit timelike vector field 
u"' (u°'Ua = —1) which is non-shearing and non- 
rotating, i.e., its covariant derivative is of the form 

V bUa ^ Ohab - UaUb, hab = Qab + UaUb, (1) 

where the acceleration and expansion 

rate 6 — VqU" are the remaining kinematic quan- 
tities of u"-; 

(ii) the Einstein tensor has the structure 

Gab = SUaUb+pgab^WUaUb+phab, (2) 
DaW = ha''VbW = 0, (3) 

i.e., the spacetime represents the gravitational field 
of either a perfect fluid with shearfree normal four- 
velocity u", pressure p -I- A and spatially homoge- 
neous energy density w — A (case S = w + p 0) 
or a vacuum {Einstein space case S — 0, where 
w = —p may be identified with A); 

(iii) the Weyl tensor Cabcd is degenerate but non-zero, 
i.e., the spacetime is algebraically special but not 
conformally flat. 

Choose null vector fields k°- and subject to the nor- 
malization condition k°'la = —1, such that 
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{qk^' + r), q>0. 



(4) 



Within the GHP formalism (cf. appendix A) based on the 
complex null tetrad m", m"), q is (-2,-2)-weighted 

and the conditions (i) and (ii) translate into 

TT + T = qK + q^^iy, X = qa, ^ -JI = q{p ~ p), (5) 
P'q-qpq^-2q{p-qp), = 5'g = (6) 



and 



$22 = — , 



$01 = *12 = $02 = 0, 

i? = 24n = w - 3p = 4u) - 35*, 
Sw = d'w = 0, P'w — qPw — 0, 



(7) 

(8) 

(9) 
(10) 



respectively. By virtue of condition (i) the magnetic 
part Hab = hvacmnC™" bdu'^u''' of the Weyl tensor wrt 
u"' vanishes [23|. In combination with condition (iii) it 
follows that the Weyl tensor is purely electric wrt 
Eab = Gacbdu'^u''' 7^ 0, the Weyl-Pctrov type is _D, and at 
each point lies in the plane S spanned by the Weyl 
PNDs (cf. appendix B for a GHP proof of these well- 
known facts). Hence, choosing fc" and 1°' along the PNDs, 
{k°-,l°-,m°-,m'') is a Weyl principal null tetrad (WPNT) 
and we have 



*o = *i = *3 = *4 = 0, 



(11) 
(12) 



Under the restrictions © and ^T^, the GHP Bianchi 
equations are given by (jA3ip - (|A36|) and their prime du- 
als. Combining these with the other equations in (|5])- P^ 
results in 



K = 


0, Z/ : 


= 0, 0- = A = 0, 


(13) 


P = 


p, p- 


= P, TT = -T, 


(14) 






P'* = -3/**, 


(15) 




= 3r*, 


g'* = -37r*, 


(16) 


P'S 


-qPS 


= S{Pq- p + qp), 


(17) 


dS- 


-rS, 


g'5 = tS, 


(18) 


P'w 


= qPw 


3S{p~qp) 
2 


(19) 



With ©-© and (Ulll-dlll) the Ricci equations, given by 
(|A25P - (|A30p and their prime duals, reduce to 



Pp = -P'p 

^, _ _ * w S* 

= -Q T + pp + TT + — + - - 

PV = V-^, 9a* = 9V = 0, 



Pp 



1. 

Tq 



, gp = g'p = 0, 



Pt = P't = 0, gr = 
—dn = 9V = 9r = — 



(20) 
(21) 
(22) 

(23) 
(24) 
(25) 



and the complex conjugates of (j24p . while the commu- 
tator relations applied to a (wp, Wg)-weighted scalar rj 
become 

[P, P']p = (wp + w,) j^rr - - + - - - j 77, (26) 
[g,g']r7=(wp-w,)('-/.p+-|-^'),7, (27) 



[P, g]?? = {-tP + pB + Wq pT)r], 
[P, d']r] = i-Tp + pd' + wp pr)77, 
[P',d]f] = {-Tp' -pd + WppT)r), 
[!>', g']?7 = i-Tp' - pd' + Wg pT)r,. 



(28) 
(29) 
(30) 
(31) 



Then the [g,g'](T), [g,g'](r), [I>,g'](T) and [P',d']{T) 
commutator relations imply 

gi/ = 2T(iJ + *-G), &H = 2t{H + -G), 

PH = p(H + F), P'H^-p{H + F), (32) 



where 



F = 2tt, G=2pp+-. 



(33) 



One checks that the integrability conditions for the sys- 
tem of partial differential equations (PDEs) are 
identically satisfied, indicating that corresponding solu- 
tions exist. Those for which is non-expanding addi- 
tionally satisfy 



M - 9P = 



(34) 
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(cf. ([Ml) and pUD)) below). However, does not follow 
as a consequence of the ansiitze; this implies the existence 
of expanding anisotropic perfect fluid models in class I (§ 
[3]). Also, the scalar invariant iip may be strictly negative, 
which is incompatible with ([M]); as a consequence, the 
class I Einstein spaces are not necessarily static f^i 12. 3p . 



2.2. Metric structure and subclassification 

The first, second and last parts of (jl3l) - P^ precisely 
account for the hypersurface-orthogonality of fc", P and 
to" o to", respectively. Thus real scalar fields u, v, (zero- 
weighted) and U, V ((—1, —1)- resp. (1, l)-weighted), and 
complex scalar fields C (zero-weighted) and Z ((1,-1)- 
weighted) exist such that 

^1/3 ^1/3 ^1/3 

daU=^^ka, daV = —^la, daC = — ^"Zq . (35) 

By (jA16p this is equivalent to 



P'u = 




Pu 


— Bu = 


B'u = 


0, 


(36) 


Pv = 




P'v 


= 9w = 


&v = 


0, 


(37) 


g'C 


= 


PC 


= P'C = 


= gc = 


0, 


(38) 




= ^'/yz, 


K 


= P'C = 


g'c = 


: 0. 


(39) 



The commutator relations (051) - (I?T|) applied to u, v, ^ 
and C then yield 



(40) 
(41) 



m = d'U = dV = d'V = 0, 
PZ ^ P'Z ^PZ = P'Z ^ 0. 



Hence, when we take these fields as coordinates, ([55)) - 
(|¥T|) imply that the zero- weighted fields UV and ZZ only 
depend on (m, u), resp. {(^, C), such that all class I metrics 
are conformally related to direct sums of metrics on two- 
spaces: 

gab = *-'/'(5ib®5S,), (42) 
g^, = 2*2/3^(^77^^^ ^2ZZ(C,C)d(„Cdh)C, (43) 
g^, = -2v|/2/3fc(J,) =_2C/F(«,«)d(„udf,)«. (44) 



The line elements of and will be denoted by ds^. 



resp. dsj 



In the case where such a two-space is not of constant 
curvature, however, we will construct more suitable coor- 
dinates in the sequel. Inspired by the GHP manipulations 
of [33*1 for type D vacua f53|i we start this construction 
by deducing suitable combinations of the scalar invari- 
ants F, G, H and From (|Jl6)) . (fTU)) and dUD-dMl) it 
is found that 

d^F^S^-i/^aa, d,G = 3*1/3^ (45) 
d,(^ = 2«'i/3xaa, d,7 = 2^-1/3^/3^, (46) 

daX = ^^/^aa, day ^^'^^ Pa, (47) 



where 

Ua =Tma + Trria, jSa = ^J.ka - pla (48) 

are invariantly-defined one-forms and 

= -^ + ^ = -^ + ^ + ^ + 2^ 

- 34>i/3 ' ^ - 3^1/3 ' ^^^1 

H + -^-G -H + 2^ + G 

' — — ■ (^0) 

Consequently, the scalar invariants 

G^3(^-x2) = 3(7-z/), (51) 
D = -x^ - Cx + F = + Cy - G. (52) 

are constant {daC = daD = 0). From and ([511) it 
follows that F, G, H and 5" are biunivocally related to 
X, y, C and I?, where 



2rr = F = x 
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(53) 
(54) 



2Aip^G-^ = y' + Gy-i?-|, 

2Wt = H = 2x^ + ix^y + Cy-D, (55) 
* = (x-|-y)VO. (56) 

Barnes [JJ partitioned class I according to the position 
of the gradient V°^' relative to S and E^. This relates 
to the vanishing of the invariants tt — —ttt or fip, maxi- 
mal symmetry of g^j^ or g^i^ and spatial rotation or boost 
isotropy of gat, as follows. 

First assume t — 0. In this case and the first 
parts of (I2S1) and (gHll-dSll) imply 

H = F = if = 0, *-G = 3a;^'2/^ 

G = -3x^, D = 2a;^ (57) 

such that X is constant. In combination with the last 
part of (HI]), (HI]) and the first parts of ([17I)-(|1H1) one gets 

rT = 0^7r = T = 0<^a; const <^ £ S. (58) 

The [9, 3'] commutator relation applied to (, ( and Z 
imply BZ ^ B'Z ^ and m' Z = 3x^^/^Z. Herewith 
the Gaussian curvature of the two-space with metric 
becomes 



{ZZy\\n{ZZ))^^-^ 
B'Z 



= _<f-2/3 



= -«'-2/3gg'(inZZ) 

-3a;, (59) 



where the dual of (|35|) was used in the calculation. In con- 
junction with the results of Goode and Wainwright [33 |. 
we conclude that ([55)) yields the class I solutions which 
are locally rotationally symmetric (LRS) of label II in the 
Stewart-Ellis classification [s^ , characterized by g^f^ hav- 
ing constant curvature = —3x. As well known (see 
e.g. the appendix of [36| ) the coordinates C and C may 
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then be adapted such that ZZ(C, C) = (1 + K^CC/'^)"'^ 
in (j43p . or an alternative form may be taken: 



1 + 



2 

^-2 



fci €{-1,0,1}. 



(60) 



Now assmne fxp =^ 0. It follows from (HOI)- (ESI), 
55| and the second parts of (|49l)-(l52|) that 



A 
3' 



(61) 



such that y is constant. In combination with psp and 
the second parts of (|Tfl) - (|48l) this implies 

^j,p = 0-^fi = p = 0-^y = const ^ G E-^. (62) 

By a similar reasoning as in the case t ~ one con- 
cludes that (|^^ yields the locally boost-isotropic Einstein 
spaces of Petrov type D, characterized by having con- 
stant curvature 



-3y, 



(63) 



such that in this case one may take UV{u,v) = (1 — 
K^uv/2y^ in (gH) and we have 

1 — 

K^ = k^Y^\ fcsG {-1,0,1}. (64) 

With and ds| written in the second form, it is clear 
that 



9./ = -*-2/3y^2 cos( Vfea;3)'d"x4 (65) 

is a HO timelike Killing vector field. 

Four subclasses of class I thus arise, which were labeled 
by Barnes as follows: 



lA : T = = ^p, IB : r = 7^ /ip, 
IC : r 7^ = pp, ID : r =^ =^ pp. 



(66) 



We proceed with the respective integrations. Notice that 
in the joint case ppr = one has 



w 



2{TT + pp)^{x + yY + K{x + yY--, (67) 

with if = for r = and K = for pp = 0. When 
r 7^ or pp 7^ we may take a;, resp. y as a coordinate, 
where (gZll-dlH]) and ^ imply 



(x + y){Tma + Trria) = d^x, 
{x + y){^ika - pla) = day. 



(68) 
(69) 



In view of and ([55]) it then remains to determine 

suitable complementary coordinates for x in g^^^ or y in 



For T 7^ 0, Frobenius' theorem and suggest to 
examine whether zero-weighted functions and / exist 
such that 



.x + y . _ _ . 

1 „ _ yjrfla - Tma) = fdaCi 
ZTT 



(70) 



By (|A16p this amounts to calculating the integrability 
conditions of the system 

Pcj) = P'ct) = 0, tS0 = -rSV = (71) 



2/ 



which turn out to be 



Pf = P'f = 0, a, W = 0. 



(72) 



These last equations have the trivial solution / = 1, for 
which a solution (f> of (|7ip is determined up to an irrele- 
vant constant. Herewith the invariantly-defined one- form 
on the left-hand side in (|70p is exact, and we take <j> as 
the coordinate complementary to x. On solving (|68p and 
(|70p with / = 1 for TOq and TOq and using (j53p we con- 
clude from (l43l) that 



dsi = — + 2Trd02 , 2tt = x^ + Cx + D (73) 

ZTT 

for classes IC and ID. Clearly, the metric solutions should 
be restricted to spacetime regions where x^ + Cx + D > 
for consistency, while 



5^" = i- 



2tt 



rd>, (74) 



x + y {x + y) 

is a HO spacelike Killing vector field (KVF). 
For pp 7^ one analogously considers 

dij} = = 0, nPij) = pP'ijj = i^-t^ (75) 

but the integrability conditions of this system are now 

gg = g'5 = 0, p,y-g^-gS^^^^. (76) 

So (7 = 1 is 072??/ a solution in the Einstein subcase = 0, 
for which we then get 



dsi 



dy^ 
2pp 



2ppdV'^ 2pp = y^ + Cy-D-^ (77) 



from dSl), (IMl), (Ei) and dZSl), with KVF 



pfc'' + pi"" 

x + y 



2pp 



(a; + y)' 



(78) 



which is timelike for pp > and spacelike for pp < 0. In 
general, the second vector field in (|78p is always HO: the 
integrability conditions of ([7S)l are checked to be identi- 
cally satisfied, such that solutions g and a corresponding 
solution ^ of (j75p exist. However, taking ■(/; as a comple- 
mentary coordinate of y eventually leads to a very com- 
plicated system of coupled partial differential equations 
for g = g(i/,ip), which is impossible to solve explicitly. 
We shall remedy this in section [3TT] but now discuss char- 
acterizing features of the Einstein space limit cases. 
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2.3. Characterizations of PE Petrov type D 
Einstein spaces 

Petrov type D Einstein spaces constitute the class 2?o 
(cf. the introduction) and are all explicitly known. The 
line elements are obtained by putting the electromagnetic 
charge parameter $o or + equal to zero in the 2?- 
metrics given by Debever et al. resp. Garcia [l5| . 

These coordinate forms generalize and streamline those 
found by Kinnersley [stI in the A = case. 

Recently, a manifestly invariant treatment of I?Oi mak- 
ing use of the GHP formalism, was presented [sll . Within 
GHP, Po-metrics are characterized by the existence of a 
complex null tetrad wrt which ([TT]) and $y — hold (i.e., 
the tetrad is a WPNT and dZ])-® with 5 = hold). Ac- 
cording to the Goldberg-Sachs theorem ([T^ holds 
and characterizes WPNTs as well. The scalar invariant 
identities (see (ssL [38l|) 

ptp = Jip, VTTf — TT, (79) 

just as ((inil-dTni), (EOll-lEII) and the first equation of ^ 
are also valid in general. From these relations it follows 
that 

ma) ^ d, (80) 

i.e., a Petrov type D Einstein space is PE if and only 
if the WPNT directions are HO. In fact, it can readily 
be shown by a more detailed analysis than in [33| that if 
the spacetime belongs to Kundt's class, i.e., if one of the 
PNDs is moreover non-diverging, one has 

/i = 0^p = 0orp-p7^07^7r + r. (81) 

Equations (4), (5) and (31) in [33] then imply 

p^-p^Q ^ p=-p^O = TT + T, (82) 

TT = —T 7^0 ^ /i = /I, p = p. (83) 

One concludes that the Kundt and Robinson- Trautman 
subclasses of Vq have empty intersection, and that the 
latter consists of PE spacetimes for which both PNDs are 
non-twisting but diverging. These results - which remain 
valid for the electrovac class V, just as the two theorems 
below - are implicit in [l4|, where the concerning PE 
metrics form the Einstein space subclasses of the classes 
labeled by 

: r = = MP, C\ : t = ^ ^^p, 
C°-: T^Q = fip, C* : T^O^pp. ^^^> 

By the Einstein space specifications S — and w = A = 
const, the boost-field q disappears from the equations 
((7))- p3p and is not determined by the geometry, in con- 
trast to the situation for perfect fiuids S ^ (cf. § 13. ip . 
Moreover, the set ([5])-®, i.e. the requirement that an 
US given by ^ exists, is decoupled from ([7))- ([55)) and is 
not needed to derive ((T5l)-(I321) from ©-([ll]). From the 
integrability of the complete set (|6))- ([33| . (|66|) and the 
above we conclude: 



Theorem 2.1 The closed set ([71)- (|33l) characterizes the 
class Vq of PE Petrov type D Einstein spaces, which 
are precisely those Einstein spaces for which the WPNT 
directions are HO or, alternatively, those which belong 
to Barnes' class I, all admitting a one-degree freedom of 
USs in all regions of spacetime. Barnes' boost-isotropic 
Kundt classes lA and IC coincide with (7*^°, resp. C°_, 
while the Robinson- Trautman members of Vq constitute 
and C* , form the Einstein space subclasses of IB, 
resp. ID, and possess non-twisting but diverging PNDs 
at each point. 

The result is in agreement with proposition IB.li which 
provides criteria for deciding when a Petrov type D 
spacetime allows for an US, regardless of the struc- 
ture of the energy-momentum tensor. The hypersurface- 
orthogonality (fT5)) - P^ of the WPNT directions corre- 
sponds to criterion 5 and is actually equivalent to a one- 
degree freedom of USs. It is worth to mention that all 
LRS II spacetimes, i.e. those exhibiting (pseudo-) spher- 
ical or plane symmetry, share this property with the Vq- 
and P-metrics. On the other hand, certainly not all PE 
spacetimes admit an US. For instance, the Godel solution 
is an LRS I PE perfect fiuid of Petrov type D, described 
in GHP by l^-l^ and 

S/2^w^p^ -3*2 = -2/ip = const. > 0, (85) 
TT = T = 0, p = qp, JI= -/i, p = -p, (86) 

q > being annihilated by all weighted GHP-derivatives; 
hence the invariant {p~Jl){p—'p) = App = Aqp^ appearing 
in criterion 2 of proposition lB.il is strictly negative, and 
it follows that the Godel solution does not admit an US. 
As another example, the spatially-homogeneous A = 
vacuum metrics 

Pi + P2 + P3 ^ pI + pj + pI = I, P1P2P3 7^ 0.(88) 

attributed to Kasner [s^ are PE [5l|. If the pi are all 
different, there is a complete group G3/ of isometries 
and the Petrov type is I. In this case dt"" is the up to 
reflection unique Weyl principal vector field and hence 
the only possible US-candidate; however, its shear tensor 
has the non-zero eigenvalues (1/3 — pi)/t and hence the 
spacetime does not admit an US. On the other hand, if 
two Pi's are equal it follows that P2 = P3 = —'2pi — 2/3 
(without loss of generality) . Then the line element repre- 
sents a Petrov type D, non-stationary, plane-symmetric 
vacuum which, according to theorem 12.11 admits a 
one-degree freedom of USs (cf. the end of this section). 

In theorem 3 of [Tl it is claimed that all vacuum space- 
times admitting an US are static, which would generalize 
Triimper's result jl^l by including Petrov type D. How- 
ever, this conclusion only holds when pp > 0. Indeed, 
a static member of Vq necessarily admits a rigid (i.e. 
non-expanding) US, such that pp > 0, cf. ([34)) . Con- 
versely, when pp = ov pp > for a PE member, it 
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admits the HO timelike KVF resp. and is thus 
static. This is in agreement with proposition IB. 31 re- 
garding fi = p = criterion 6" tells that in fact all boost- 
isotropic spacetimes, with tt — —t wrt a WPNT, are 
static, while for fip > one checks that criterion 2" is sat- 
isfied by virtue of ([TT]) - ([55)) . In appendix B the freedom 
of the rigid USs and HO timelike KVF directions (static 
observers) in these cases is also specified, which is in ac- 
cordance with a result by Wahlquist and Estabrook [i^l ■ 
In summary we have: 

Theorem 2.2 A Petrov type D Einstein space is static 
if it admits a rigid US. This is precisely the case when 
the spacetime is PE and has a positive or zero scalar in- 
variant p,p, being the product of the divergences of (non- 
twisting) Weyl principal null vectors and l'^ subject 
to k'^la — —1. For pp > there is an up to reflection 
unique rigid US, deflned from the geometry by Q and 
q ~ fi/p and parallel to the unique HO timelike KVF di- 
rection. For pp = 0=>p = p — (classes lA and IC) all 
USs are rigid USs and have a one-degree freedom, while 
the HO timelike KVF directions are parametrized by two 
constants. 

For completeness we display standard coordinate forms 
of the PE Petrov type D Einstein space metrics, as re- 
covered here by (gH), ^ and (Uni), (El, CI), dZB- 

C°° corresponds to and From dST]), dSH), 

(pTjl and one deduces that 

K-^ = -3x ^ -3y = , ^2 = = 4:X^ 0. 

Rescaling (, u and w by a factor (2a;) one arrives at 

, 2 2dCdC 2duAv A - „ . , „ 

1 -I- f CC 1 - fww 2 

This represents the Einstein space limit <f>o of 
Bertotti's static and homogeneous electrovac family with 
cosmological constant (4ll.l4^. exhibiting spatial rotation 
and boost isotropy (complete group Gq of isometrics). 
The A = limit yields fiat Minkowski spacetime. 

C°+ and C°_ correspond to and ([771), resp. ([M]) 
and ([731). Making use of replacing in the C°+ (C°_) 
case the coordinate y [x) by r = — (2m)^/'^/ {x-\-y), rescal- 
ing the remaining coordinates by a factor {2m)~^/^ and 
writing Y = (2m) ^/"^ > one finds 

ds^ = fd^2 ^ ScosiVkifdTf) + - 5gk{r)Ax^, 

^ ' 9k(r) 

Off! A 

gk{r) = k -r\ k = X(2m)i/3 e {-1, 0, 1}, 

r 3 

with (5 = 1 for C"+ (^ = -1 for C"_). These solutions 
have a complete group G4 of isometries acting on space- 
like (timelike) three-dimensional orbits, and for A = 
correspond to Kinnersley's case I (IV) with I = 0. The 
static region of (gfc(r)>0) yields class A in the clas- 
sification of static Petrov type D vacua by Ehlers and 



Kundt i; C0_ is static everywhere and corresponds to 
class B. Regarding C^, the subcase k = 1 reproduces 
after ^ i-^ 7r/2 — ^ the well-known forms of the spher- 
ically symmetric Schwarzschild-Kottler interior and ex- 
terior metrics fH, the subcase A; = A = 0, r>0 
{9k{^) < 0) gives another form of the plane-symmetric 
Kasner metrics (cf. supra). 

C* corresponds to ((75)) and (1771) . which gives the line 
element 



1 



dx^ 



{x + yY \f{x) 



d?/2 



9{y) 



A 



f{x)=x' +CX + D>Q, g(y) = -/(-y)--. (89) 

The KVFs d^"" and 9^" generate the complete, abelian 
group G2 of isometries. For A = 0, ([55]) is the form 
of the C-metric obtained by Levi-Civita and recovered 
by Ehlers and Kundt, and corresponds to Kinnersley's 
case III A. It is generally assumed - and suggested in the 
original paper [6| - that the Kinnersley- Walker form 



1 



hm4>^ + 



k{ri)d'ip' 



h{^) = 1 - ^2 2ma^^ > 0, k{r]) = -h{-r]) (90) 

equivalently describes the gravitational field of the A = 
C-metric. However, this is not entirely correct. Equating 
the Lorentz invariants appearing in the right hand sides 
of the equations in ([5n | - ((5^ . calculated for the metrics 
([89| and yields 



X = ~{2my^ ( Q!^+ ) , y = -(2m) 3 (a?] - 

Dm / V Dm 



G = -- 



1 



3(2m)^ 



1 



54m2 



(91) 



Hence only covers the range C < 0, £> > -2 (-y) ^ , 
whereas in general the constant scalar invariants C and 
D are allowed to take any real value. Yet, the cubic ](x) 
has discriminant — 4C^ — 27D^] thus it has three distinct 
real roots if and only if 



C < -3 



^ C<0, |D|<2 



(92) 



Thus (|?T|) is compatible for this case, and by further 
rescaling (j) and ijj with a factor a(2m)~^/'^ one arrives 
at ([90]) ; (|92p is equivalent with ma < 1/ -\/27, leading to 
the physical interpretation of two uniformly accelerating 
masses. Recently, Hong and Teo [45| introduced a nor- 
malized factored form for this situation, which greatly 
simplifies certain analyses of the C-metric. A further 
coordinate transformation can be made such that the 
Schwarzschild metric is comprised as the subcase a = Q. 
This was further exploited for the full D-class in [l9| . 

Finally, we write down the equations which determine 
all USs for a member of C*, in the coordinates y and 
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of 



Let 



V5(y) 



in the static region and 



U^^-{x + y)./^)dy-, V^^ 



x + y 



in the non-static region, and gauge-fix = ([/" -|- 
V'')/yf{2), l" = {U" - V'')/^{i). The unit timehke 
field @ is an US if and only © holds; this translates to 

1 = q{y,i-') and 

g{v){q ± l)q^y + {qT + g' {v)q) = 0- (93) 

Here and below the upper (lower) signs should be taken 
in the static (non-static) region. For solutions q — q{y)^ 
i.e. q^^ = 0, direct integration of ([M)) yields 



9{y){q{y)Tlf ^E±q{y), 



(94) 



with i?+ > and < constants of integration. No- 
tice that in the static region the solution q{y) = 1 yields 
the unique static observer. In the case q^^ 7^ the so- 
lutions get implicitly determined by an equation of the 
form ip = ip(y,q)^ on applying the method of character- 
istics for first-order PDEs (see e.g. H^). In the subcase 
where C — D = A = 0, g{y) reduces to y^ and this 
equation reads 



^ = - 



y(gT 1)^ 



3q2 



(95) 

Here Z is a free function of its argument, making the 
one-degree freedom of USs more explicit. Replacing ([75)) 
by ([SO)) does not alter these equations, i.e., the above 
remains valid for C°+. Then C = £> = is equivalent 
to a; = if-L = 0, cf. dST]) and ([59]); for A = the non- 
static region y < corresponds to the plane-symmetric 
Kasner vacuum metrics, where y = —(3^/2)"^/'^ and a 
rescaling of the other coordinates recovers (|57)) . the USs 
being determined by (|94)) and (j95|) with the lower sign. 



3. PERFECT FLUID GENERALIZATIONS OF 
THE C-METRIC 

3.1. Line element 

We resume the integration of class I started at the end 
of § 12.31 We thereby focus on the subclass ID charac- 
terized by r 7^ 7^ [xp. Let us first summarize what we 
did so far. We started oflF with the closed set 
of first-order GHP equations in the seven (weighted) real 
variables 5*, S*, w, /i, p, g, 3't and the complex variable r. 
These variables are equivalent to two dimensionless spin 
and boost gauge fields, e.g. r/r and [ij p. and seven real 



scalar invariants. The boost and spin gauge fields could 
serve to invariantly fix the tetrad - the ID members be- 
ing therefore anisotropic - but can be further ignored. 
For the C*-Einstein spaces, S — ^ and w = A = const, 
and we remarked that q is not a part of the intrinsic de- 
scribing set of variables. Hence we end up with four real 
scalar invariants in this subcase. These invariants are 
equivalent to the two constants C and D and two inde- 
pendent functions x and y, which we took as coordinates 
and in terms of which, on adding two coordinates (p and ^ 
related to the symmetries, the corresponding C*-metric 
can be expressed. In the perfect fluid case S 0, ^ 
gives the boost field q which, starting from an arbitrary 
gauge (fc", 1°"), turns given by Q into the invariantly- 
defined fiuid four- velocity. The four invariants and their 
use persist, just as the coordinate (f>. However, the scalar 
invariants w and S are no longer constant and ^ is no 
longer a suitable coordinate. Thus we need one more 
scalar invariant for our description and one remaining 
coordinate complementary to y. 

For the first purpose it is natural to look at the kine- 
matics of the fluid, which are fully determined by 



Here 



(gr - /") 



(96) 
(97) 

(98) 



is the intrinsic spacelike vector field, which determines at 
each point the up to reflection unique normalized vector 
orthogonal to and lying in the PND plane S, while 
Mil and are the component along w", resp. projection 
onto of the acceleration u". In analogy with ((96|) we 
define the invariant 



b = 2V(a^;c)"^°m.'^. 
The relation with GHP quantities is 



2q 



(99) 



(100) 



iH\^i2qr'/'{P'q + qPq) = ^-b, (101) 



-U^ = TTUa + Tnia ^ Ota = 



Notice that pOOp is equivalent to 



hUa- bVa = pka - pL = Pa 



Pq_ 

x + y' 

_ day 

x + y' 



(102) 



(103) 



In combination with (US)) -(EH), (fTO^ and (fTUH)) imply 



2tt = iitu^'" = + Cx + D, 



2pp -b^ = y^ + Cy- D 



w 



(104) 
(105) 
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We choose b as the final describing invariant and use b 
and U|| as auxihary variables. In view of (|100p - (|102p one 
deduces that the difi'erential information for S, w and b 
comprised in is precisely 



DaS = 


-SUa, 




(106) 


daW = 


-u{w)Ua, U(w) 


~3bS, 


(107) 










dab = - 


-u{b)ua, u{b) = 


-v(b) + - 


-^)-f,(108) 


v(6) = 


-^(3y^ + C). 




(109) 



From (|109p it follows that b is non-constant, such that 
we may see the second part of (jlOSp as a definition oi u\\. 
P07p is nothing but the energy resp. momentum conser- 
vation equations for a perfect fluid subject to DaW — 0. 
The first part of equation (|108l) confirms that Da9 ~ 
P, [l^l , whilst the second implies again that the expan- 
sion scalar does not vanish in general (cf. the end of § 12.11 
and below). 

For the second purpose we rely on the hypersurface- 
orthogonality of u"^ by assumption: zero-weighted real 
scalar fields t and / exist such that 

dj - lua. (110) 
The integrability condition hereof is 

Dal = -IUa = -I{Ua +U\\Va), (111) 

which is equivalent to 

m^Tl, v(/) = (112) 



From 5^0, pU5)) and plU)) it follows that t is function- 
ally independent of y (and of x and cf)) and we take it as 
the fourth coordinate. With the aid of ([TTUl)- pTT1) . pIIS)) 
and the first parts of p07p and pOSp precisely tell that 
A = w and b only depend on t. Hence b — b{y,t) 
from (fTU5)) . On using (fT02 | - (fT03)) . the first part of (fTT2)) 
is equivalent to J = J{y,t), where J = ^yj^- Eliminat- 
ing -uii between the second parts of (|108l) and (jll2p . and 
using v{x + y) = —b{x + y) imphed by (|102p - (|103p . yields 

b^^^{J) ^bJu{b)+A{x + y). (113) 
Inverting and pTU]) we get 

{x + y)ua = bJdat, [x + y)va = bJdat (114) 





or dually 



+ 69/, =59/. (115) 

x + y bJ x + y y 



Thus in the chosen coordinates (jll3p reads 

dyJ{y,t)^b{y,tr^b'{t)-A{t)]. (116) 



From (gH), dZSl), g^, = i^aVb - UaUb)/{x + yf and the 
only remaining equation (|107p we obtain the line element 



ds^ = {x + y)-^[dsl + dsl], 

dsi = + 2Trd0^ , 2rr = x^ 
Itt 

ds\ = (bJdt - -{ bJdt", ^ 



where 

b^b{t), w^w{t), A 



bJS 



(117) 
Cx + D, (118) 

(119) 



A{t), (120) 



2{x + y) 

w'{t) = Qb{t)A{t) ^ daW = 66 Adat = ?,bSua, (121) 
b = b{y, t) = ^/y^ + Cy-D + 6(t)2 - w{t)/3, (122) 



J = J(y,t)-[fe'(i)-A(i)] 



dy 



+ L{t), (123) 



with L{t) a free function of integration. The solutions 
are defined and regular in the coordinate regions 

2tt = x^ + Cx + D > 0, (124) 
b{y, tf = y^ + Cy - D + b{tY -w{t)lZ>Q. (125) 

Notice that we nowhere used S ^ Q explicitly in the 
above integration procedure. Therefore, the above line 
element describes the complete class ID, including the 
C*-vacuum limits which correspond to w{t) — A and 
A{t) = 0, cf. (|120p . In this case the coordinate transfor- 
mation {t,y,x,(l)) ^ (V':2/j2^j0)j which connects (jll7p - 
p23p to the original form eliminates b{t) and L{t) 
and follows from Q, (fTUO)) and ^T^, giving 



'^o.W = (Aifca + Plaj 



2^p 

Jdat + Z-7Z 



(bUa - bVa^ 



day- (126) 



Hence, ip — ip(y, t) and it is the solution of the consistent 
system 



dtijj = J, dyip = —r: 



b b^-b^ 



(127) 



the integrability condition hereof being precisely (jll6p 
with A{t) = 0. The transformation is singular at degen- 
erate roots of &2 and at the union of the black hole and 
acceleration horizons [ia,il] ^2 _ 52 = f(^_y-^j_/^/^ ^ g, 
which separate the static from the non-static regions. Let 
us emphasize that the &(<)-freedom is essentially a free- 
dom in the choice of coordinates. The form ((89)) de- 
scribes the full C-metric manifold; y can take any value, 
and the sign of f{—y) -|- ^ is positive in the static re- 
gion and negative in the non-static region. In the form 



(|117p - (lll9p y is always spacelike and we have constructed 
t as a synchronized timelike coordinate corresponding to 
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an US with associated expansion rate 9{t) = 3b{t); 
for fixed b{t) the range of y is constrained by (jl25p and 
only this subregion of the manifold is described by the 
coordinates. E.g. pT7)l - p^ with A{t) ^ 0, w{t) = A/3, 
b{t) — and L{t) — 1 (which formally reduces to ([55)) 
on putting t — ip additionally) only describes the static 
part of the C-metric, the vector field then lying along 
the unique HO timelike KVF direction. However, in 
the neighborhood of any point with coordinate label y, 
the metric can be described by pi7l) - (jl23l) . by choosing 

K0'>/(-y) + A/3. 

We neither used r ^ 0. This implies that the line 
element of the complete class IB, characterized by ^ 
= T and constituted by all LRS H Einstein spaces and 
shear-free perfect fluids with DaW = 0, is described by 
PT7| - (fn5)) . with (fTT5| replaced b y CT . This class was 
first described by Kustaanheimo [i^ and rediscovered 
by Barnes [1], both using different coordinates (see also 
(16.49), (16.51) in [16]). 

Of course, the result (|117p - p23p could have been ob- 
tained without referring to GHP calculus. Barnes [i| 
showed that the metric can be written in the form 



As' 



with f ^ f{x). Indeed, from ([TTl) . ([TTC)l and ([T^ it 
follows that {x + y)va is exact: (x + y)va = daz; z is 
used as a coordinate instead of y, and one puts Jb = , 
Z = Z{y,t). Notice from pM)l that now 



y^Y{z,t), 



e = 3Yte- 



(129) 



Let us directly attack the field equations in these coor- 
dinates, thereby correcting [l|. One can check that only 
four of the field equations are not identically satisfied (the 
indices 1 to 4 label the Weyl principal tetrad vectors nat- 
urally associated with (|128D ) : 

G3i = 0= -Ytz + YtZ,,, (130) 
Gn - G33 = = (131) 

/' - 2r,, + ix + Y) {zl + Z - /72) , 
Gii=p= (132) 

2 {Y,u - Y^tZ.t) - Y,,Z,, - f 12) (x + Y) 

+ (z 3, + z^^) {x + r) V 3 (y2 + / _ r2g-2z^ ^ 
G33 G44 = s = (133) 

2{x + Y) [Y,, - Y,Z^, + (Yu - YtZ^t)] ■ 

Hence, if supplemented with 9 ^ Y_t ~ 0, these equations 
are the ones obtained in Barnes |il.]: equation (|130p = 
v{9) — was missed out, and both equations (|132l) and 
([155)1 differ from equations (4.23), resp. (4.24) in [Ij] by 
a term 2{x + Y)Y^tte~''^ ■ Thus, it is clear that with 
these differences a correct non-expanding solution can 
be found, but the analysis of expanding solutions will be 
incorrect. 

Differentiating ()132p twice wrt x yields d'^J{x)/dx'^ = 
0, whence 



Substituting this in equation ()132p . and equating coeffi- 
cients of powers of x, leads to 

Z ^^(z, t) + Z^^{z, tf = iaY{z, t) - b, (135) 
F,,(z,t) - ^-Y{z,t)b+^aY{z,tf. (136) 

The solutions Y[z,t) of the last equation are defined by 

^Y(z.t) ^ 

-z + f2{t) = 0, (137) 



■y/ ar^ — br^ + cr + fi (t) 



which can be solved for z in terms of Y . This eventually 
suggests to transform coordinates from {z,t) into {y,t), 
with y — Y(z,t). Rescaling and translating coordinates 
allows us to set a = 1 and b — 0. One can check that the 
remaining equations lead exactly to equations (|116p and 
dnj), recovering solution ([TT7| - ([n5)) . 



3.2. Properties 



Consider the metric ()117l) - ()123p . for which we assume 
henceforth that it describes a perfect fluid {A{t) 7^ 0). In 
contrast to the Einstein subcase, is now the unique 
invariantly-defined fluid velocity, and the expansion rate 
9{t) and energy density w{t) — A of the fluid are scalar 
invariants. Expressions for the pressure p -I- A and the 
components and U|| of the acceleration follow from 

(Una, (unHD-cMi), mB and ^im -. 



ibJ)iy,t) 



x + y f3y^ + C b'{t)-A{t) 



b{y,t) 



ibJ)iy,t) 



The fluid is non-shearing and non-rotating, i.e. w° is an 
US. Because of ([T5)) - ([H)) criterion 5 of proposition IB.ll 
is satisfled, such that there is a one-degree freedom of 
USs. These can be found by taking g = 1 in ([4]) and 
([M)l . hereby flxing the (fc", Z°)-gauge geometrically, and 
solving (O with q replaced by Q, the USs then being 
{Qk'' + P)/v/2Q. This yields Q = Q{y,t) and 



bJ[b{Q + 1) + b{Q - l)]Q,y + {Q~ l)Q,t 
= -2Q{Q-l)b{bJ),y 
Q{Q-l 



f{x) = ax^ + bx^ + cx + d. 



(134) 



{3y' + C)bJ + 2{b' - A) 



If the class is to be used as a cosmological model, 
it is interesting to discuss the intrinsic freedom. By 
(fT^ and ([TTi)) we have that 2A{t)dat = Sua and 
J{y,t)dat — {x + y)ua/b are invariantly-defined one- 
forms, and hence so is L{t)dat because of (|123p . It follows 
that ■j{t) is a scalar invariant. Moreover, as A(t)dat is 
exact we may remove the only remaining coordinate free- 
dom on t by putting A{t) = 1, such that the conservation 
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of energy equation (|121l) can be considered as a defini- 
tion 9{t) = w'{t)/2. Hence, in this most general pic- 
ture for S* 7^ 0, the scalar constants C , D and invariants 
^{t), w{t) characterize the model within the class. Notice 
that the presence of two invariantly-defined, distinguish- 
ing free functions could have been predicted, since after 
elimination of , there are two scalar invariants u(&) and 
u{S) remaining unprescribed in the system of equations 
(fTU5)i - (fM)D . 

In this fashion however, the physical implications re- 
main obscure: it would be nice to have a free func- 
tion, with a clear physical interpretation, instead of 
L/A. Spacetimes with L{t) = have w{t) as the only 
free function. If L{t) ^ 0, L(t) can alternatively be 
fixed to 1 by a t-coordinate transformation. In this 
case the metric structure functions display the expan- 
sion scalar, the energy density and the pressure (since 
1A{t)Aat = (w + p)ua); these are related by energy con- 
servation (|12ip . where w{t) and A{t) can be chosen freely. 
Alternatively, one can subdivide further in = and 
6 0. In the case 9 = 0, the energy density w — A 
is constant because of (I12ip and can be chosen freely, 
just as A{t). In the most interesting case 9 0, w{t) 
and 9{t) can be chosen freely, determining Sua via p2ip . 
Thus class ID provides a class of anisotropic cosmologi- 
cal models with arbitrary evolution of matter density and 
(non-zero) expansion. 

Regarding symmetry, all perfect fluid ID models admit 
at least one KVF d^'^ given by ([7^ . which at each point 
yields an invariantly-defined spacelike vector orthogonal 
to u-*"" and lying in S"*- . If </> is chosen to be a periodic co- 
ordinate, with range given by [—ttE^ttEI, the spacetime 
is cyclically symmetric. We will then refer to the region 
F = f{x) = 0, where the norm of 9^" vanishes, as the 
axis of symmetry [Z^l [E^. Finding the complete group 
of isometrics and their nature is trivial in our approach. 
The functions x, y, w and L/A are invariant scalars, such 
that K^daX = K'^daV = K'^daW = K'^da^ = for any 
KVF if. As the ID models are anisotropic, it follows 
that the complete isometry group is at most G2 , and if it 
is G2, both w and L/A are constant. Conversely, when 
w and L/A are constant we have 9 = 3b = from (|12ip . 
b = b{y) from (fT^ and J{y,t) = -A{t)F2(y) from ([^5)1 . 
By redefining the time coordinate such that A{t) = 1 one 
sees from PT71) - P^ that dt" is a HO timehke KVF. We 
conclude that the ID perfect fluid models have at least 
one spacelike KVF 8^°" , which may be interpreted as the 
generator of cyclic symmetry. They admit a second in- 
dependent KVF if and only if both scalar invariants w 
and L/A are constant, in which case the spacetimes are 
static and the complete group of isometrics is abelian G2, 
generated by 8^°" and dt°' . 

Consider the case where f{x) has 3 real non-degenerate 
roots Xi, i.e. ((5^ holds. If xi < X2 < X3 then f{x) > 
for all X S ]a;i,a;2[- Furthermore, we let </) be a periodic 
coordinate. The ratio between circumference and radius 
of a small circle around the axis, a; = xi or a: = a;2, is 



given by 



2TTEy/f{^) 



= -ttE {3x1 + C) 



(138) 



respectively 



lim 



27r£;V7M 



= ttE (3x1 + C) 



(139) 



It is only possible to choose the parameter E such that 
the complete axis is regular, if 3xf + C = ^ (3a;2 + C"). 
However, eliminating C and D between this equation and 
f{xi) = f{x2) = implies xi = X2. Consequently, if 
f{x) has three real non-degenerate roots, the spacetime 
contains a conical singularity. This echoes the properties 
of the C-metric |43|, and suggests the presence of a 
cosmic string. 



4. CONCLUSIONS AND DISCUSSION 

A new class of Petrov type D exact solutions of Ein- 
stein's field equation in a perfect fluid with spatially ho- 
mogeneous energy density has been presented. It con- 
sists of all anisotropic such fluids with shearfree normal 
four- velocity, and generalizes a previously found class to 
include non-zero expansion. The analysis and integra- 
tion was rooted in the 2+2 structure of the metric and 
use of invariant quantities. This approach clarified the 
link with the vacuum C-metric limit, and certain prop- 
erties of the vacuum case are inherited. However, the 
presence of the perfect fluid defines generically two extra 
invariants. For the expanding solutions, this translates 
into an evolution of energy density and expansion which 
can be chosen freely. This subclass contains only one 
(potentially cyclic) symmetry. 

The viability of these solutions as a low-symmetry 
class of cosmological models is subject to further re- 
search. More in particular, it should be clarified whether 
a thermodynamic interpretation of the perfect fluid can 
be made [50j - it is certainly not possible to prescribe a 
barotropic equation of state p = p{w). The relation with 
the C-metric also suggests to further examine the arising 
coordinate ranges, properties of horizons, and whether 
an interpretation as a perturbation for small masses of a 
known PF solution exists for certain members. 
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Appendix A: Geroch-Held-Penrose (GHP) 
formalism 

The GHP formalism (TgI . \2T\ is a complex, scalar for- 
malism, which is a 'weighted' version of the Newman- 
Penrose (NP) tetrad formalism. Use is made of a complex 
null tetrad (ei°, 62°, 63°, 64°) = (to°, to°, Z", fc"), where 



fc"L = -1, m°-m„ = 1 



(Al) 



and all other inner products vanish. To put it in other 
words, at each point one takes a timelike plane, two vec- 
tors fc" and lying along its real null directions, and 
two vectors and to" lying along the complex con- 
jugate null directions of the orthogonal spacelike plane, 
these pairs of vectors satisfying the normalization con- 
ditions (jAip . We use the labels a, b etc. for the tetrad 
indices. The basic variables of the formalism are the spin 
coefficients (P-g- = ea°Vc(eg)Q = -Pg^g) 

K = P414, T = P413, 0- = P4ii, p = P412, (A2) 
f^ = r233, 7r = P234, A = P232, M = r231, (A3) 

the 9 independent components of the traceless part of the 
Ricci tensor Sab = Rab - jRgab, 

*oo = ISabk^'k', $22 = ^Sabn", (A4) 
$01 = ISabk^'m", $12 - ^Sabl^'m", (A5) 
$02 = ^Sabrn'^m'', $n - ^Sab {k^l' + m^m") , (A6) 



with — (f>io , the multiple 



n 



24 



(A7) 



of the Ricci scalar, and the 10 independent components 
of the Weyl tensor Cabal, 

■^Q^Cabcdk^m^k^m'', ■^i^Cabcdl^mH'^m^ (A8) 
*i - Cabcdk''l''k'm^, *3 = Cabcdl^kH^m'', (A9) 
*2 - Cabcdk^rn'm^l''. (AlO) 

Changes of the tetrad leaving the null directions 
spanned by l'^, m" and m° invariant, and at the same 
time preserving the normalization conditions (lAip . con- 
sist of boosts 



and spatial rotations 

Quantities transforming under (jAll|) - (jA12p as 



(All) 
(A12) 

(A13) 



are called well- weighted of type {wp,Wq) or (wp,Wq)- 
weighted {zero-weighted in the case of type (0,0)). They 
have boost- weight wb{ti) — and spin- weight 

{t]) — . One can check that the GHP basic vari- 

ables are well-weighted, their weights following from the 
definitions (IAToI) and (|ATT|) - (|AT3| - see also equa- 
tion (7.36) in [l||. E.g. wb{i') = -2, ws(r/) = -1, im- 
plying ly is of type (-3,-1). The following derivative oper- 
ators are defined such that a well-weighted quantity is 
transformed in a well- weighted quantity: 

Dal] = ea(r/) + P34aWB(?7) V + Pl2aW5(ry) 77. (A14) 

When T] is of type (wp, Wg) one can check that 
^BiDaV) = WB(??)+WB(a), wsiDAv) = ^s{ri)+wsia), 
where 

{1, a = 4, r 1, a = 1 

-1, a = 3, , ^(a) = < -1, a = 2 
0, a = 1,2 [ 0, a = 3,4. 

One uses the notation 

d = Di, d' = D2, P'^Ds, P = Di. (A15) 

Notice that the differential of zero- weighted scalars / can 
be expressed as 

daf = -P'fka-Pfla + &fma+Qfrna (A16) 

= -l(/)fc,-k(/)L + m(/)TO„ + m(/)m, (A17) 

= -u{f)Ua + ^{f)Va + m(/)TO„ + m(/)TOa, (A18) 

where and w° are related to k°- and according to 
dH), resp. IMl). 

The basic (or 'structure') equations of the GHP for- 
malism are 

(a) the commutator relations of the weighted derivatives, 
in the joint Da notation given by 



[Da, V = 2Pf..ji?g,; + ws(r;) [R^,,~, + '2T ,,[aT\,^^) V 

+Wsiri) (^Rj^2ab + 2ric[arj=2|gj) V + W B{b)T j,4,dDyr] 

-\-wsib)Ti2aD'^r] - WBia)T^^i^Da7] - ■ws{a)T^^i^Dari; 

(b) 12 complex Ricci identities (or 'equations'), namely 

with [ab] — [14], [23] (the complex conjugates corre- 
sponding to [ab] = [24], [13]); 

(c) the Bianchi identities (or 'equations') 



^[f(^\db\cd\> 



'"^^abelc^y] + ^l[cRdf]eb H[c^df]ea- 



One can show that, after writing the directional deriva- 
tives ea in terms of the weighted derivatives Da, these 
basic equations (a)-(c) form a consistent, closed system 
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of PDEs in the variables (|X2|) - (IA10|) and with formal 
derivative operators Da- Compared to the NP formalism, 
the 6 complex Ricci identities which concern directional 
derivatives of the non- well- weighted NP spin coefficients 
a, /3, 7 and e (corresponding to [ah\ — [12], [34]) have 
been absorbed in the commutator relations (1). Explic- 
itly, for a (wp, w^)-weighted scalar one gets 

[i>,p'](77) = (^ + T)g(77) + (^ + T)g'(,7) 

■^{kV - TTT + n - $11 - 1'2)Wp 77 

-f(7«7- TTf + n - $11 - *2)w,r?, (A19) 
[6, g'](ry) = (M- 71)1^(7?) + (p-p)I>' (77) 

+ (A(T - Tip - n - $11 + *2)Wp ?7 

-(A^-/ip-n-$ii +l'2)w,ry, (A20) 
[P,g](ry) = 7fP(ry) - nP' {r]) +-pd{ri) + a& {r^) 
+ {k^ - an - 'I'i)wp 77 

-h(KA-7fp- $oi)wg?7, (A21) 

together with the equations obtained by applying the 
complex conjugate and/or prime dual operation to (jA21l) . 
This prime dual operation is generated by interchanging 
k"' •(-)■ /" and O m", which comes down to 



K -n- —V, T O — TT, fT o —A, p —pL, (A22) 
$2-j2-j, 0*4-4, (A23) 

Pop', go g'. (A24) 



The interchange (jA24|) means that (Prj)' — P'r/' etc., and 
is due to (IA14I) and 



-wb(?7), ws(?7') = -ws(?7), 
i.e. Wp(?7') = -Wp(77), W5(r/) = -Wq{7]). 

Regarding complex conjugation one has Prj = Prj, dr/ = 
d'rj and 

wb{v) = wb(77), W5(77) -ws{t]), 



I.e. w. 



,(77) = W,(77), Wq(77) = Wp(77). 



Explicitly, the 12 complex Ricci identities read 

Pt - P'k = (t + n)p + (t + TT)a + $01 + (A25) 

gp - gV = (p - p)t + - + $01 - , (A26) 

Pcr-gK= (p + p)cr + (7r-r)K + «'o, (A27) 

Pp- g'K = p^ + (TCT - KT + KTT + $00, (A28) 

PV - gr = -cr/i - Ap - + Ki7- $02, (A29) 

P'p - gV = -pp - Xa-TT + Kv - 2n - I'a (A30) 



and their prime duals (IA25P '- (IA30P '. Finally, the Bianchi 
identities involve weighted derivatives of the Riemann 
tensor components. In full generality they are given in 
ref. [Hi, (7.32a-k), or [H, (4.12.36-41). 

The formalism is especially suited for situations where 
two null directions are singled out by the geometry, such 
that fc" and 1°" can be chosen along them. In particu- 
lar, the Weyl tensor of a Petrov type D spacetime has 



precisely two PNDs; choosing and 1°' along them is 
equivalent to condition (jlip . and a complex null tetrad 
realizing this condition is called a Weyl principal null 
tetrad (WPNT). When (0) and ^ are both satisfied, 
the Bianchi identities reduce to 

= ct(2$ii -f-3^'2) A$oo, (A31) 
P*2 + I^'^oo + 2Pn = p(2$ii + 3«'2) - p^oo, (A32) 
P$ii + P'$oo + 3Pn = 2(p + p)$ii - (p + p)$oo,(A33) 
^2 + 2gn = -t(2$ii - 3*2) + j7*oo, (A34) 
g$ii - 3gn = 2(r - 7f)$ll - 77$oo + ^$22, (A35) 
g$oo = k(2$ii - 3^-2) - 7f$oo (A36) 



and their prime duals (|X3T|) '- (|M6l) '. 

In general, the GHP formalism may be used to find a 
class of solutions, defined by a particular set of properties. 
One first translates these properties in terms of GHP 
variables, yielding (algebraic or differential) constraints 
on the system of basic equations, then recloses the result- 
ing extended system (integrability analysis) and finally 
describes the corresponding metrics in terms of coordi- 
nates (integration). These coordinates are four suitable, 
functionally independent zero- weighted scalars /; they 
may be combinations of (derivatives of) basic variables, 
appearing in the reclosed system S itself, or 'external' co- 
ordinates associated to HO vector fields due to Frobenius' 
theorem. The geometric duals of the null tetrad vectors, 
and hence the metrics gab — —'2k(ah) + '^i^iaJ^b): &re 
obtained by inverting (jA16|) for the chosen /'s. Eventu- 
ally the remaining equations of S are written in terms 
of these coordinates and the resulting PDE's are solved 
as far as possible. We refer to [5l| for enlightening dis- 
cussions, and to e.g. [52| or this work for illustrations. 
In particular for Petrov type D spacetimes, notice that 
zero-weighted combinations of WPNT spin coefficients 
and their weighted derivatives (e.g. pp or d'r) are scalar 
(Lorentz) invariants x, which are thus annihilated by any 
present KVF K", K(x) = Ckx = 0. This facilitates 
the detection of KVFs. More generally, zero- weighted 
tensor fields Tab..., algebraically constructed from the 
Riemann tensor, WPNT vectors and covariant deriva- 
tives thereof, are invariantly- defined by the geometry, 
and C-KTab... = 0. 



Appendix B: (Rigid) shearfree normality and 
staticity of Petrov type D spacetimes 

Consider (an open region of) a spacetime and a unit 
timelike vector field u° defined on it. Choose a null vec- 
tor field k°-. At each point, fc" and span a timelike 
plane E, the first null direction of which is spanned by 
fc". Construct the null vector field 1°' by taking at each 
point the unique vector lying along the second null direc- 
tion and satisfying k°'la = — 1. Then is decomposed 
as in dH), where q ^ , A ^ -{V^k^Ua)'^ ■ The field 
defined in ([M)) determines at each point the up to 
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reflection unique unit spacelike vector lying in E and or- 
thogonal to u'^. The electric and magnetic parts of the 
Weyl tensor wrt can be decomposed as 



Eab = CacbdU^U = (^'2 + *2)[WaW6 - m(^amb)] 



^'4 + (7^*0 ^*3-g*l 



2q 



-marrib 



2q 



l{aVb) 



c.c, (Bl) 



Hab = ^^^C"'^tdu'u'' - - 'i'2)[VaVb " m(,mb)] 



*4 



2 



2q 



rriamb + 2 j= — m(^a'Vb) 



2q 



c.c.(B2) 



If u° exists such that Hab = 0, the Weyl tensor is purely 
electric (PE) wrt u", the spacetime itself being also called 
PE. A criterion in terms of Weyl tensor concomitants, 
deciding whether this is the case, follows from the flow 
diagram 9.1 in ^16y] and theorem 1 in (53| . 

Suppose now that the spacetime admits a unit timelike 
vector field satisfying ([T]) - corresponding to an US, i.e. 
forming the tangent field of a shearfree and vorticity-free 
cloud of test particles. Within the GHP formalism based 
on k"' and 1°" as introduced above, this is the case if and 
only if a (-2,-2)-weighted field q exists satisfying ([S])-®. 
By virtue of these relations, the [3, 5'] (5) commutator 
relation yields (HH), adding 2q[ (|M6)) - ((X2B)) '[_to the [P' - 
qP,d']{q) commutator relation gives "ifs + q'^i = and 
the combination (f^^ - fM?)) ' + q[(M2il' - (P^ ] 
produces - q^^^o = 0. Hence Hab = from (IB2]) 



'4 — q — 'J- nence Hab 
and if we choose fc" to be a (multiple) PND, = 
(^Tq = \iJ^=0), then also 1°- is a (multiple) PND, ^-4 = 
{^4 = V^a =0). Hence the spacetime must be either 
conformally flat (all zero), and then USs are always 
admitted (see e.g. (6.15)) in [16]), or purely electric (PE) 
and of Petrov type D or /, the Weyl tensor being PE wrt 
For Petrov type /, there are 4 distinct PNDs, and u'^ 
is the up to reflection unique timelike vector lying along 
the intersection of the planes spanned by two particular 
pairs of PNDs. For Petrov type D, k°- and I"' can be 
taken to be the multiple PNDs, and lies in the plane 
S spanned by them. 

Propositions 4 in (s^ and 16 in [59] imply intrinsic, 
easily testable criteria for deciding when a Petrov type I 
spacetime admits an US, resp. is static. Here we present 
likewise criteria in the Petrov type D case. These crite- 
ria are invariant statements, in terms of GHP basic vari- 
ables and weighted derivatives associated to an arbitrary 
WPNT. Given a Petrov type D spacetime in coordinates, 
the determination of the PNDs, and hence the WPNTs, 
is straightforward and can be performed covariantly. It 
then suffices to fix one WPNT and calculate the appear- 
ing spin-boost covariant expressions by using definitions 
(IM)) - (|XTO)) and (|XT4)) . For complex (2fc, 2fc) -weighted 
scalars (fc e Z) z — Re(z) + i Im(z) we mean with z > 
(z < 0) that z is real and strict positive (negative) in the 
sequel. 

It turns out that, given (fTT|) - P^ . the integrability con- 
ditions of ^ are identically satisfied. Thus we find: 



Proposition B.l A Petrov type D spacetime admits an 
US if and only if, wrt an arbitrary WPNT, ^2 is real and 
one of the following sets of conditions holds: 

1. fj 7^ 0, the scalar invariant Act > 0, and go = ^1'^ 
satisfies ([S])-®; 

2. p ^'p^ the real scalar invariant {p, — Jl){p — 'p) > 0, 
and qo = —{p — 'p)/{p — p) satisfies ©-([S]); 

3. X — (7 — p — ]I = p — ]} — 0, the scalar invariant 
Kiy and one of the following situations occurs, 
where go defined in each subcase satisfies ([SJ and 
where b = (tt + t)/k, c = v/k: 

3a. Im(6)Im(c) > and qo = Im(c)/Im(6) also 
satisfies — Re(6)go -|- Re(c) = 0; 



3b. b ^ Re(fe), c<Oandqa = {b + -4c)/2; 

3c. & > 0, c > 0, &^ > 4c, and qo = {b + 
Vb^ - Ac)/ 2 or qo={b- Vb^ ~ 4c) /2; 

4. X — a — p — 'p = p — 'p = 0, and either 

4.1 K — ^ {n + t)v > 0, and go = I'/iir + t) 
satisfies or 

4.2 K =^ = ly, (tt + t)k > and go = (tt + t)/k 
satisfies ©; 

5. the WPNT directions are HO, i.e., ([I3])-([I11) holds. 

The subdivision of case 3 stems from a straightforward 
analysis of the second equation of ([5]). In cases 1, 2, 
3a, 3b and 4 there is a unique US, whereas there may 
be one or two USs in case 3c. Due to the number and 
nature of the equations © there is a one-degree freedom 
of USs in case 5, where the condition — ^'2 can 
be dropped since it is implied by the imaginary part 
(IA301) ' and 



of (|A30P + (|A30P ' and ([131) -((Til). Important examples 
of spacetimes satisfying criterion 5 are the Petrov type 
D purely electric Einstein spaces and their 'electrovac' 
generalizations (see [38, 56] and § 12. 3p and all spacetimes 
with (pseudo-) spherical or planar symmetry (which 
constitute the LRS II Lorentzian spaces, see [33)- These 
examples all satisfy ([7]) on top of ([T^ - (fH)) and are 
further characterized by $00 — $22 = ($11 =) 0, resp. 
Tr = T = m = (cf. [34]). 

The spacetime will admit a unit timelike vector field 
u° satisfying 



Ua-b 



-UaUb, 



(B3) 



corresponding to a rigid US or modeling a rigid non- 
rotating cloud of test particles, if and only if a (-2,-2)- 
weighted field q exists satisfying (O-® and ((34)) . Notice 
that, given ([M)) . the third equation of ([5]) is identically 
satisfied. Hence we have 

Proposition B.2 A Petrov type D spacetime admits a 
rigid US if and only if, wrt an arbitrary WPNT, ^2 is 
real and one of the following sets of conditions holds: 
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1'. condition 1 with the third equation of ^ replaced 

by (121; 

2'. the scalar invariant /jp > and qq = fi/j) satisfies 

3 '-5'. conditions 3-5 with fj, — JI — p ~J) — replaced by 
p = p = 0. 

In case 5', the spacetime possesses geodesic, shearfree and 
non-diverging PNDs {K = a = p — 0, i' = X = fj, = 0)- 
thus belonging to Kundt's class - and HO Weyl principal 
complex null directions {X — a — t: + t — 0), and admits 
a one-degree freedom of rigid USs. 

The spacetime is static if and only if it admits a 
HO timelike KVF. An equivalent characterization was 
given by Ehlers and Kundt |5|: the spacetime is static 
if and only if a unit timelike vector field exists for 
which shear, vorticity and expansion scalar vanish, i.e. 
(|B3p holds, and for which the acceleration m° is Fermi- 
propagated along the integral curves of u": 

U[aUb] = 0. (B4) 

The field m° is then parallel to a (HO and timelike) KVF 
and identified with a congruence of static observers. By a 
long but straightforward calculation, thereby simplifying 
expressions by means of ([S])-®, dMl), (!A25|) . (|A25P ' and 
the [P, !>']((/) commutator relation, one shows that the 
extra condition (jB4[) is equivalent to 

{qK + q'^V){Pq + y^) - 2PV + 2qpT 

+$12 - q^Qi = 0, (B5) 
PPq = TTT + TTT — q{KTr + 'kW) — q~^{yTX + Fvr) 

+2$u-^ + 2^'2. (B6) 

In case 5' above, the Ricci equations (jA25|) . (jA28p and 
(E28)) ' yield Pt = $oi and $oo = $22 = 0, and so (|B5)) - 
(IB6p reduces to 

$i2+'Z$oi=0, (B7) 

PI>q = -2Tr + 2$n-4 + 2*2. (B8) 



In the subcase $01 = $12 = of (lB7|, the \P,P'\, [I>,S] 
and [P, 5'] commutators applied to q yield 

P'Pq= -qPPq+i:Pqf, (B9) 
5f-g = Tpq, g'Pg = Tpq. (BIO) 

The compatibility requirement of (jB8p - (|B10p with the 
commutator relations for Pq gives the single condition 

P'R + qPR^^. (Bll) 

According to the Sach's star dual of the LRS criterion 
in [Hi, the subcase P' R = Pi? = of (|BTT]) precisely 
corresponds to a boost-isotropic spacetime with vr + t = 
0. From the above we conclude: 

Proposition B.3 A Petrov type D spacetime is static if 
and only if, wrt an arbitrary WPNT, one of the following 
sets of conditions holds: 

l"-4". ^'2 is real, conditions l'-4' hold and go additionally 
satisfies ((B5|) - (|B6)) : 

5"a. condition 5' holds, the scalar invariant $oi'l'2i < 
and go = -$i2/$oi satisfies dH) and ((B8| : 

5"b. condition 5' holds, (f>oi = $21 = 0, the scalar in- 
variant (P'i?)(Pi?) < and go ^ -P' RjPR satisfies 
© and (iBSl) : 

6". the spacetime is (locally) boost-isotropic and tt -|- 
r = 0. 

The HO timelike KVF directions are parametrized by 
two constants in case 6" [60], are 1 or 2 in number in 
case 3"c and are unique in all other cases. 



[1] A. Barnes, Gen. Rel Grav. 4, 105 (1973) 

[2] A. Weyl, Ann. Phys. (Germany) 54, 117 (1917) 

[3] T. Levi-Civita, Rend. R. Acad. Lincei, CI. Sci. Fis. Mat. 

Nat. 26, 519 (1917) and 27, 183, 343 (1918) 
[4] I. Robinson and A. Trautman, Proc. Roy. Soc. Land. A 

265, 463 (1962) 
[5] J. Ehlers and W. Kundt, in Gravitation: an introduction 

to current research, edited by L. Witten, page 49 (New 

York: Wiley, 1962) 
[6] W. Kinnersley and M. Walker M, Phys. Rev. D 2, 1359 

(1970) 

[7] J. Bleak and B. Schmidt, Phys. Rev. D 40, 1827 (1989) 



[8] V. Pravda and A. Pravdova, Gzech. J. Phys. 50, 333 
(2000) 

[9] W. B. Bonnor, Gen. Rel. Grav. 15, 535 (1983) 
[10] K. Hong and E. Teo, Glass. Quantum Grav. 22, 109 
(2005) 

[11] J. B. Griffiths and J. Podolsky, Glass. Quantum Grav. 

22, 3467 (2005) 
[12] H. Farhoosh and R. L. Zimmerman, Phys. Rev. D 21, 

2064 (1980) 

[13[ J. Bicak and V. Pravda, Phys. Rev. D 60, 044004 (1999) 
[14| R. Debever, N. Kamran R. G. McLenaghan, J. Math. 
Phys. 25, 1955 (1984) 



16 



D. A. Garcia, J. Math. Phys. 25, 1951 (1984) 
H. Stephani, D. Kramer, M. A. H. MacCallum, C. 
Hoenselaers and E. Herlt, Exact Solutions to Einstein's 
Field Equations (Second Edition) (Cambridge: Cam- 
bridge University Press, 2003) 

J. F. Plebariski and M. Demianski, Ann.Phys. (N. Y.) 98, 
98 (1976) 

E. T. Newman, L. A. Tamburino and T. Unti, J. Math. 
Phys. 4, 915 (1963) 

J. B. Griffiths and J. Podolsky, Exact space-times in Ein- 
stein's general relativity (Cambridge: Cambridge Univer- 
sity Press, 2009) 

J. N. Goldberg and R. K. Sachs, Acta Phys. Polon., 
Suppl. 22, 13 (1962) 

J. J. Ferrando, J. A. Morales and M. Portilla, Phys. Rev. 
D 46, 578 (1992) 

J. J. Ferrando, J. A. Morales and M. Portilla, Phys. Rev. 
D 50, 2567 (1994) 

M. Triimper, J. Math. Phys. 6, 584 (1965) 
M. Triimper, Z. Phys. 168, 55 (1962) 
A. Barnes, J. Phys. A 5, 374 (1972) 

H. Stephani, Comm. Math. Phys. 4, 137 (1967) and 5, 
337 (1967) 

R. Geroch, A. Held and R. Penrose, J. Math. Phys. 14, 
874 (1973) 

J. Podolsky and J. B. Griffiths, Phys. Rev. D 63, 024006 
(2000) 

J. Podolsky, Czech. J. Phys. 52, 1 (2002) 
O. J. C. Dias and J. P. S. Lemos, Phys. Rev. D 67, 
064001, 084018 (2003) 

R. Emparan, G. T. Horowitz and R. C. Myers, JHEP 
0001 007, 021 (2000). 

R. Penrose and W. Rindler, Spinors and Spacetime vol 1 
(Cambridge: Cambridge University Press, 1986) 
S. B. Edgar, A. Garcia-Parrado, J. M. Martin-Garcia, 
Class. Quantum Grav. 26, 105022 (2009) 
S.W. Goode and J. Wainwright, Cen. Rel. Grav. 18, 315 
(1986) 

J. M. Stewart and G. F. R. Ellis, J. Math. Phys. 9, 1072 
(1968) 

P. Szekeres, Comm. Math. Phys. 41, 55 (1975) 
W. Kinnersley, J. Math. Phys. 10, 1195 (1969) 
S. R. Czapor and R. G. McLenaghan, J. Math. Phys. 23, 
2159 (1982) 

[39] E. Kasner, Amer. J. Math. 43, 217 (1921) and Trans. 
A.M.S. 27, 155 (1925) 



[40] H.D. Wahlquist and F. B. Estabrook, J. Math. Phys. 7, 
894 (1966) 

[41[ B. Bertotti, Phys. Rev. B 133, 1331 (1959) 

[42[ I. Robinson, Bull. Acad. Pol. Sci. Ser. Sci. Math. Astron. 

Phys. 7, 351 (1959) 
[43[ K. Schwarzschild, Preuss. Akad. Wiss. Berlin, Sitzber., 

189 (1916) 

[44[ F. Kottler, Ann. Phys. (Germany) 56, 410 (1918) 
[45| K. Hong and E. Teo, Glass. Quantum Grav. 20, 3269 
(2003) 

[46[ Y. Pinchover and J. Rubinstein, An introduction to par- 
tial differential equations (Cambridge: Cambridge Uni- 
versity Press, 2005) 

[47[ J. B. Griffiths, P. Krtous and J. Podolsky, Class. Quan- 
tum Grav. 23, 6745 (2006) 

[48[ P. Kustaanheimo, Comment. Phys. Math., Helsingf. 13, 
8 (1947) 

[49[ H. Bondi, Mon. Not. Roy. Astron. Soc. 142, 333 (1969) 
[50| B. Coll and J. J. Ferrando, J. Math. Phys. 30, 2918 
(1989) 

[51[ S. B. Edgar, Gen. Rel. Grav. 12, 347 (1980) and 24, 1267 
(1992) 

[52[ S. B. Edgar and G. Ludwig, Gen. Rel. Grav. 29, 1309 
(1997) 

[53[ C. B. G. Mcintosh, R. Arianrhod, S. T. Wade and C. 

Hoenselaers, Class. Quantum Grav. 11, 1555 (1994) 
[54[ J. J. Ferrando and J. A. Saez, Class. Quantum Grav. 14, 

129 (1997) 

[55[ J. J. Ferrando, J. A. Morales and J. A. Saez, Class. Quan- 
tum Grav. 18, 4939 (2001) 

[56[ R. Debever and R. G. McLenaghan, J. Math. Phys. 22, 
1711 (1981) 

[57[ More precisely equations (48)- (53) of ref. [33l [. There are 
some typos in these equations: the integers 3 and 9 
should be omitted in (50) and (51)-(52), respectively, 
while there should be 9k instead of k in (53). 

[58[ If one of the pi's is zero then so is a second one; this 
reproduces Minkowski. 

[59[ Strictly speaking, with the terminology of [l^ . the space- 
time has in general only a cyclic symmetry, as the axis 
will be shown to be irregular, and consequently not part 
of the spacetime. 

[60[ The two constants result from the integrations of (|B8[1 - 
(|B10|) . giving Pq, and consecutively of 



